Pulse-induced acoustoelectric vibrations in surface-gated GaAs-based quantum devices 
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We present the results of a numerical investigation which show the excitation of acoustoelectric 
modes of vibration in GaAs-based heterostructures due to sharp nano-second electric-field pulses 
applied across surface gates. In particular, we show that the pulses applied in quantum information 
processing applications are capable of exciting acoustoelectric modes of vibration including surface 
acoustic modes which propagate for distances greater than conventional device dimensions. We 
show that the pulse-induced acoustoelectric vibrations are capable of inducing significant undesired 
perturbations to the evolution of quantum systems. 

PACS numbers: 85.35.Gv, 73.23.-b, 77.65. Dq 



INTRODUCTION 

Quantum computers have the potential to efficiently 
solve "hard problems" , for which there are no known ef- 
ficient classical algorithms jjy . The realization of a quan- 
tum processor is currently a major challenge at the fore- 
front of experimental physics. Semiconductor quantum 
dots form the basis of many solid-state devices being de- 
veloped for the implementation of quantum logic gates. 
They are believed to be promising candidates for pro- 
ducing scalable systems, as they take advantage of the 
technology available for the fabrication of semiconductor 
technology. The essential requirements for the realiza- 
tion of such devices are the coherent manipulation of an 
electronic degree of freedom in each dot and the control 
of coupling between adjacent dots. Since the theoretical 
work of Loss and DiVincenzo [2] , experimental progress 
has been made_in the realization of spin qubits ^ and 
charge qubits [4] based on quantum dots formed by the 
surface gating of a two-dimensional electron gas (2DEG) 
where voltage pulses are used to control exchange interac- 
tion and tunnelling, respectively, between adjacent dots. 
However, the scaling of such systems requires the identi- 
fication and minimization of the decoherence processes, 
which arise from the coupling of the quantum dot to its 
environment. This requires a precise understanding of 
the dynamical evolution of both the relevant electronic 
degree of freedom and the environment to which it is 
coupled. 

In our recent theoretical work 0, we had imple- 
mented mechanical boundary conditions in the form of 
time-dependent traction forces on the surface of a GaAs- 
based device, to excite acoustoelectric vibrations. By 
utilizing the piezoelectric nature of GaAs, time- varying 
electric fields form the basis of surface acoustic wave gen- 
eration through interdigitated transducers. The time- 
varying electric fields such as those used to manipulate 
spin and charge qubits in quantum information process- 



ing applications may also excite acoustoelectric vibra- 
tions. These vibrations would be undesirable for the op- 
eration of sensitive quantum dot systems. 

In this paper, we show through the numerical solu- 
tion of the equations of motion in GaAs, the excitation 
of acoustoelectric modes of vibration in GaAs-based het- 
erostructure devices, due to pulsed electric fields. We 
discuss the results in relation to the construction of scal- 
able quantum dot devices. 



THEORY 

The equations of motion in a general piezoelectric ma- 
terial are Isl 
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where Ui is a displacement vector, (j) is the electrostatic 
potential field, c^j.^ are components of the elastic tensor, 
efj are components of the material dielectric (the super- 
scripts E and S indicate measurements under constant 
strain and electric field respectively), g is the mass den- 
sity, and Cijk are components of the piezoelectric tensor. 
The mechanical boundary conditions require the pre- 



scription of the traction forces Ci 



at each external 



surface, where ri is a unit vector perpendicular to an 
external surface, and cr.y are components of the stress 
tensor, defined by 

(Jij = -SkijEk + cfjkieki, (3) 

^here E, = and e., = i(|f- + The electro- 

static boundary conditions require the prescription of the 
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FIG. 1: (a) A schematic representation of the device used 
in our simulations. The orientation of the crystal axes in 
relation to the coordinate a:xes is shown. The parameters a 
and d represent the width of the gates and distance between 
the gates, respectively, (b) The shape of the pulse used in the 
electrical pulse experiments. The parameters A, tr and td are 
the amplitude, rise time and pulse duration respectively. 



normal component of the electric displacement at the free 
surface of the medium. The electric displacement vector 
is defined by 



D, 



^ij^j + eijk£jk- 



(4) 



The equations consist of a system of coupled linear par- 
tial differential equations second order in space and time. 
The solution of these equations for non-trivial boundary 
conditions requires numerical methods. Our numerical 
solution method, which utilizes a finite element approx- 
imation for the spatial part and a second order finite 
difference discretisation in time, is suitable for the sim- 
ulation of electromechanical waves We have used 
eight-noded brick elements corresponding to linear basis 
functions Ni^ resulting in an overall spatial and temporal 
convergence rate of 2 for the error. An operator splitting 
strategy is employed to split the coupled Ui-tj) problem . 
The numerical formulation and verification is described 
in more detail in Ref. [^. 

For GaAs, the non-zero components of the piezoelec- 
tric tensor have the value ei4 = —0.16 Cm^^. The 
non- vanishing components of the permittivity tensor are 
en = 622 = £33 = es = 13.18eo (eo = 8.85 • IQ-^^ Fm"!), 
and the non-vanishing components of the elastic tensor 
are c^^^^ = Cyyyy = Czzzz = cu = 11.88 • 10^° Nm"2, 
Cxxyy = Cyyzz = Czzxx = ci2 = 5.38 • 10^° Nm"^, and 
Cxyxy = Cyzyz = Czxzx = C44 = 5.94 • 10^" Nm"2. All 
other non-zero components of the elasticity tensor can be 
determined by applying its symmetry properties Cijki — 
Cjiki — Cijik = Ckiij ■ The mass density of GaAs g has the 
value of 5.36 • 10^ Kgrn'^ 0. 

Figure [TJa) shows a schematic representation of the 
system we have modelled in our simulations. The co- 
ordinate X and z -axes are aligned with the crystal [Oil] 
and [100] directions, respectively, as is frequently the case 



in the realization of GaAs-based quantum dot systems. 
For this orientation, the properties of acoustoelectric vi- 
brations are well documented 0, [lH and hence it is also 
convenient for the interpretation of our results. 

The anisotropic nature of GaAs will result in some de- 
pendence of the acoustic properties on the crystal ori- 
entation. However, this is expected to be small as the 
components of the piezoelectric tensor are equal and the 
components of the elasticity tensor are the same order 
of magnitude. In general, quantum dot applications will 
include electric fields across components in the orienta- 
tion investigated in this paper. The presence of a thin 
layer of AlGaAs is neglected for simplicity, and has been 
shown to have little effect on the acoustic properties 0]. 
As in our previous work we avoid a Hartree-Fock type 
self-consistent calculation as the benefits of this are likely 
to be small compared to drawbacks associated with vast 
computational resources. 

In our simulations, we implement a trapezoidal form 
for the applied voltage pulse as shown in Fig. ^h). 
Trapezoidal pulses are implemented in experiments 
where coherent manipulation or measurement of single 
electron charge and spin -systems have been achieved 
[E S, 0] ■ Figure Wljo) shows the geometrical features of 
the pulse. In experimental demonstrations, the param- 
eters of the pulse are such that tr << 0.1 ns, 0.1 ns < 
td < 1000.0 ns. We apply the voltage bias to the central 
gate relative to the outer gates. The gate width and sep- 
aration are chosen to be those frequently used in experi- 
ments to define quantum dots; a = 250 nm d = 250 nm. 
The mechanical gates are excluded from our simulations 
as they are capable of perturbing the induced vibrations 
. The presence of mechanical gates has been shown to 
reduce the amplitude of surface acoustic modes but also 
affect its vibrational frequency. To model the applied 
voltage bias, we apply Dirichlet boundary conditions on 
the regions of the surface which the gates would occupy 
and a Neumann-type condition elsewhere through setting 
Dz = 12]. For the mechanical boundary conditions, we 
apply traction-free boundary conditions i.e. cr.y ■ h ~ 
on the external surfaces. 

The precise characterization of the dependence of the 
effects of a pulse on all the parameters of this system is 
beyond the scope or intention of this paper. Our aim is to 
show the effects on the underlying lattice of electric-field 
pulses frequently used in quantum dot experiments. We 
first demonstrate the effects of a "typical" pulse used in 
experiments, and then consider the dependence of these 
effects on the parameters of the applied pulse, t^ and tr- 



RESULTS AND DISCUSSIONS 

Figure ^a,) shows the oscillations of the electric po- 
tential (f), along the x axis and 100 nm below the surface 
(as this is the standard depth of a 2DEG) , resulting from 
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a pulse with A = 1 V, tr = 0.025 ns and td = 0.3 ns. 
From time t = 0.4 ns, the waves propagate outward from 
the center at velocities ranging from ~ 2750 ms^^ to ^ 
5000 ms~^, with the lower velocities likely to correspond 
to surface acoustic waves and the higher velocities to 
other modes which include bulk waves. 




FIG. 2: (a) Time series plot of the electric potential as a 
function of x coordinate from time t — 0.4 ns to t = 2.0 ns, 
with At = 0.1 ns. The induced waves propagate outward from 
the center at a; = 8 /im with increasing time. For clarity, the 
inset shows the curve at t = 2.0 ns with a maximum amplitude 
of ~ 2.0 mV. The data was taken at a depth of 100 nm from 
the surface, (b) Grayscale plot of the pulse induced electric 
potential at time t = 2.0 ns. The regions of higher electric 
field at (a; ~ 3 /im) and {x ~ 13 fim outlined) correspond to 
single minima (or maxima) of a surface acoustic wave. 

As the inset of Fig. [HJa) shows, the maximum ampli- 
tude of the waves at the time t = 2.0 ns is ~ 2.0 mV. 
The linearity of the underlying equations of motion would 
require a linear relationship between the amplitude of 
the applied pulse, and that of the pulse induced oscil- 
lations. The pulse induced vibrations are small com- 
pared to the direct electrostatic coupling as well as rel- 
evant electrical characteristics of conventional quantum 
dot devices including the charging energy and energy- 
level spacings. In particular, the pulse-induced vibrations 
would be screened in the conductive regions. However, 
in the unscreened or depleted regions, the pulse induced 
vibrations will be visible. Of particular relevance, is in 
applications where electrical pulses are implemented to 
align electronic energy levels between a lead and a sin- 
gle quantum dot or between adjacent quantum dots. In 
such systems, the pulse-induced vibrations could momen- 



tarily perturb the alignment and lead to uncertainty in 
the evolution in the electronic system. 

Moreover, surface acoustic waves with amplitudes less 
than 0.01 mV have proved to provide an efScient means 
of implementing quantized charge pumping [l3| . 

The acoustic vibrations could be described as phonons 
which have been shown in experiment [l^ and theory 15 1 
to couple strongly to electrons. When the energy level 
separations in a quantum dot are equal to the energies 
of the excited phonons, the phonons could be absorbed 
by electrons in the dot. In a double quantum dot sys- 
tem, this could lead to inelastic tunnelling between the 
adjacent dots [l^ . 

Figure m^a) also shows the ability of the surface acous- 
tic wave modes to propagate for at least four microns 
along the surface with no decay in amplitude. The men- 
tioned four microns is only a limit of the computational 
domain and we would expect the induced surface acous- 
tic wave to propagate over much longer distances. Fig- 
ure IHb) shows grayscale plots of the induced electrical 
potential at time t = 2.0 ns. The stronger potential fields 
at (a; 3 /im) and (a; ^ 13 /im) are clearly localized at 
the surface (i.e. z = 4 /im) with the amplitude at two 
wavelengths into the depth (< 10 /iV) substantially re- 
duced compared to that at the surface (^ 1.0 mV), as is 
characteristic of surface acoustic waves in GaAs. 

In the experimental realization of quantum informa- 
tion processing devices [1, 01 the applied voltage pulses 
are of the order of 5 mV. Our simulations predict the 
excitation of surface acoustic waves with amplitudes of 

10 /iV, which is comparable to the inelastic tunnelling 
energies and the energy level bias applied to a double dot 
system [i, 17|. To demonstrate an effect of this pulsed- 
induced surface acoustic wave on a adjacent double quan- 
tum dot system, we solve numerically the relevant time- 
dependent Schrodinger equation; 



(5) 



where e{t) is the energy difference between the uncoupled 
localized states, A{t) is the energy difference between the 
stationary states when the local charge states are degen- 
erate, and ai are the Pauli matrices. In our calculation, 
e{t) is computed by assuming that the quantum dot cen- 
ters of the adjacent system are separated by a distance 
of 200 nm and have a zero energy-level bias in the 
absence of acoustoelectric vibrations. This double-dot is 
located 5 /im laterally away from the source of surface 
acoustic waves. We compare the state of the neighbor- 
ing double-dot system with IvPp) and without I'i'o) the 
electric potential wave, to determine the influence of the 
acoustic vibrations. Figure[3]shows the infidelity between 
the states, as a function of A chosen to correspond to that 
in experiments (and is assumed constant for simplicity) 

:4]. 
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FIG. 3: The effect of a pulse- induced surface acoustic wave 
on a double quantum dot device. 



The oscillations of Fig. [3] approximately commensurate 
with the frequency of the acoustic vibrations (~ 3 GHz) . 
The infidelity 1 - K^pl^fo)!^ osciUates between and 0.6 
which corresponds to a substantial perturbation to the 
evolution of the quantum state. This result is simply 
representative of what may happen, showing that signif- 
icant errors can occur for typical parameters and that a 
judicious choice of parameters can mitigate the error. 

In a system consisting of an array of quantum dots on 
the surface, each of which implements electric-field pulses 
to control an electronic degree of freedom, we would ex- 
pect the excitation of surface acoustic waves, which would 
be capable of traversing through multiple quantum dots. 
The necessity to incorporate the information regarding 
the electric- field pulses applied to every quantum dot, 
into the description of the dynamics of each quantum dot, 
potentially poses a significant challenge for the scalabil- 
ity of these systems for quantum information processing 
applications. 

In addition, the readout of the final state of a quan- 
tum system requires the preparation, manipulation and 
measurement of an ensemble of identically prepared sys- 
tems and hence, the repeated application of an electrical 
pulse. As a result, we could expect a substantial accumu- 
lation of acoustoelectric vibrations within the underlying 
lattice. 

The dependence of the pulse-induced acoustoelectric 
vibrations, on the duration of the pulse td is investigated 
by varying td from ns to 1.0 ns, with tr — 0.025 ns and 
A — 1.0 mV. The measurements correspond to 0.5 ns 
after the end of the pulse train. Figure HJa) shows 
the computed root mean square (rms) and maximum 
of the modulus (max{modulus)) of the induced oscilla- 
tions computed at 100 nm below the surface for data 
between x — nm and x — 1& /ini, for each value of 
td- Both measures are included as the rms takes into ac- 
count all wave modes (surface and bulk) whereas from the 
max{modulus) we see the maximum amplitude i.e. only 
the surface mode. One can see a rapid increase of these 
measures up to td ~ 0.18 ns, after which the rms becomes 
flat at ~ 0.3 mV, while the max{modulus) measure sta- 
bilizes at 1.2 mV. The dramatic change in behavior at 
td '~ 0.18 ns, can be understood by examining the defor- 



mation of the underlying lattice. Figure [H^b) shows the 
z component Uz of the lattice displacement as the volt- 
age pulse is applied. The x component exhibits similar 
behavior is not therefore shown. As an external electric- 
field pulse is applied, the lattice slowly deforms. After a 
certain amount of time tc (with the voltage bias on) , the 
lattice deformation reaches its equilibrium value (for the 
given electric potential), and then oscillates about this 
position before stabilizing. The time tc is ~ 0.18 ns for 
the device modelled in this work. We would expect the 
time tc to depend on the material, the device geometry, 
crystal orientation with respect to the coordinate axes 
and the model for the applied electric field (which in our 
case avoids a full self-consistent calculation). 

In our simulations, the response time tc of the lattice 
to the applied electric fields is much greater than the 
frequently applied rise times tr of the applied pulse which 
in the ideal case is zero. Hence, for the regime of interest 
in semiconductor quantum dot applications, variation of 
the rise time of the pulse would have little effect on the 
pulsed- induced vibrations. 

We may also conclude from our investigations of vary- 
ing the rise time and duration of the pulse, that the pulses 
implemented in quantum information processing appli- 
cations excite acoustoelectric vibrations of the maximum 
amplitude (for a given amplitude for the applied pulse). 

Future work may involve the investigation of novel 
schemes which may reduce the effect of the acoustoelec- 
tric vibrations. For the present time, we point out that 
surface acoustic waves may be utilized for quantum in- 
formation processing applications [l^. Moreover, mate- 
rials which have a spatial inversion center such as Si are 
immune to the excitation of acoustoelectric vibrations 
through electric field pulses and thus have a significant 
advantage for quantum control and quantum computa- 
tion. Devices based on such materials are currently being 
developed for quantum computation 3, 2^, l2l| . 

Many spin-readout proposals rely on mapping a spin 
state to a charge state, since the detectors that are likely 
to be used are Quantum Point Contacts (QPC's). There- 
fore. the results are also relevant to spin qubit systems 



CONCLUSION 

In conclusion, we have demonstrated the excitation of 
acoustoelectric modes of vibration in a GaAs lattice as 
a result of sharp nano-second electric-field pulses applied 
to surface gates. The excited wave modes include surface 
acoustic modes which have been shown to travel at least 
four microns along the free surface of a material. We have 
demonstrated the substantial perturbation these waves 
may have on the evolution of quantum systems. Whilst 
the effect we have considered does not decohere the sys- 
tem in the normal sense of becoming correlated with the 
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FIG. 4: (a) The dependence of the amphtude of the pulse- 
induced acoustoelectric vibrations <j), on the pulse duration td- 

(b) Time series plot of the z component of displacement 
across the device for a sequence of times t = 0tot = 0.5ns 
with At = 0.025 ns. A small offset is applied between curves. 

(c) For clarity, the value of «z at a; = 8 for the sequence 
of times t = to t = 0.5 ns. 

environment, the fact the acoustoelectric waves induced 
by gates at one place can significantly affect the state 
of a remote double-dot system, means that the results 
are very significant for quantum information processing 
applications. 
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